Energy gaps of Hamiltonians from graph Laplacians 
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The Cheeger inequalities give an upper and lower bound on the spectral gap of discrete Laplacians 
defined on a graph in terms of the geometric characteristics of the graph. We generalise this approach 
and we employ it to determine if a given discrete Hamiltonian with non-positive elements is gapped 
or not in the thermodynamic limit. First, we define the graph that corresponds to such a generic 
Hamiltonian. Then we present a suitable generalisation of the Cheeger inequalities that overcomes 
scaling deficiencies of the original version. By employing simple examples we illustrate how the 
generalised Cheeger inequalities can successfully identify gapped or gapless phases and we comment 
on the computational complexity of this approach. 
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Introduction: In physics, the variational approach 
plays a central role in theoretical and numerical approx- 
imations of spectral properties of Hamiltonians. When 
applied to the first excited state it can give an upper 
bound to the spectral gap of the Hamiltonian, provided 
that the ground state is known. Unfortunately, a lower 
bound to this gap is not readily available. Such a bound 
could help in understanding fundamental properties of 
many body systems that depend on the behaviour of the 
spectral gap in the thermodynamic limit [THH]- A direct 
approach to this question is, in general, a hard computa- 
tional problem as it requires the diagonalisation of expo- 
nentially large matrices. It is, thus, fascinating that the 
spectral gap of discrete Laplacians defined on a graph can 
be upper and lower bounded by general geometric char- 
acteristics of the graph [U [5] . This is achieved with the 
help of the Cheeger constant that determines the extent 
a bottleneck configuration appears in a graph. 

Here we derive a new version of the Cheeger inequali- 
ties that can estimate the energy gap of stoquastic Hamil- 
tonian, i.e. Hamiltonians with all off-diagonal elements 
real and non-positive. For that we need to first associate 
a graph Laplacian to the Hamiltonian provided we know 
its ground state. This is not a major drawback as there 
are large families of physically relevant states, e.g. the 
Matrix Product States, that are ground states of Hamil- 
tonians which are not known to be gapped or not in two 
or higher dimensions [BJ. Such an important example is 
the two-dimensional AKLT model that can support uni- 
versal quantum computation by measurements only, but 
it is not proven yet if it is gapped which would establish 
its fault-tolerance. 

The upper Cheeger bound is based on the familiar to 
physicists variational method. To generalise the lower 
bound to the case of Hamiltonians we are forced to de- 
viate from the usual methodology [5J and employ the 
duality between the maximum flow and the minimum 
cut of a graph [7]. This provides much greater flexibil- 
ity and allows us to efficiently lower bound the spectral 
gap of a Hamiltonian in the thermodynamic limit. We 



also employ an additional structure, the reduced graph, 
that allows to optimise the lower bound beyond what is 
currently possible [5]. 

The methodology presented here sets the framework 
to study the spectral gap of Hamiltonians in terms of 
geometric characteristics of graphs without any approx- 
imations. Hence, the difficulty in estimating the energy 
gap of a Hamiltonian, when its ground state is known, 
has been translated to finding the Cheeger constant of 
the corresponding graph, which is known to be an NP- 
complctc problem [8]. Our approach aims to give a new 
and general methodology that can facilitate the evalua- 
tion of energy gaps theoretically or numerically [3 . 

Hamiltonians and Graph Laplacians: We shall 
consider Hamiltonians, H, that are Hermitian N x N 
matrices with real non-positive entries. By employing 
the Perron-Frobenious Theorem [9] it can be shown that 
the smallest eigenvalue Ao of H is negative and all other 
eigenvalues are strictly larger in absolute value. Hence, 
we can denote the N eigenvalues as 

Ao < Ax < A2... < Ajv_i. 

Moreover, the eigenvector ip — (cvi)^g 1 corresponding 
to A is unique and all its components, ojj, are positive 
satisfying ^ \cti\ 2 = 1. We are interested in estimating 
the spectral gap of H given by the difference between the 
lowest and second lowest eigenvalues 

A(ff) = Ai-A . (1) 

The first excited state tpi — {fli) 1 ^ 1 has components 
that are both positive and negative in order for the or- 
thogonality condition, (^qjV'i) = 0j to hold. 

From the Hamiltonian, H, we next define the Lapla- 
cian operator L 

L = -X l N + D^HD=h 0+H - i= j: (2) 

where D is the diagonal matrix D — diag(ao, ...,Q!jv-i) 
and 11 n is the TV-dimensional identity. In contrast to 
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the usual paradigm 5J Laplacian L is a non-symmetric 
matrix that has a number of useful properties. Its 
rows sum up to zero Lij = — Ao + J2j Hij a j/ a i = 
— A + XQOti/oti = 0] and its lowest eigenvalue is with 
left eigenvector it = (ajj, a 2 N _ l ), i,e, 7ri = 0, and 
right eigenvector D _1 ?/>o = (1, l) 7 ^ = 1, i.e. LI = 0. 
Most importantly L has the same energy gap as H, 
i.e. A(L) = Ai - Ao = A(£T) [L(U _1 |Vo}) = and 
Lp- 1 |V'i» = (Ai-A )(^- 1 |V'i»]- 

We now show how to construct a graph from the 
Laplacian L ^ [TO]. Consider a set of vertices V = 
{0, N — 1}. For each pair of vertices which do not 
need to be distinct, the edge has a positive weight 



(3) 



Hence, the set of edges E of the graph involve all 
for which ^ 0. As we have determined the vertices 
and the edges of the graph we can now define the degree 
of each vertex in the following way 



h E H -. 

3&V 



Q; 



|A |af. 



(4) 



Thus, we could think of |Ao| as the 'bare' degree of i and 
di as the 'dressed' degree of i. Similarly, —Hij could be 
thought of as 'bare' weight of the edge (i, j) and Wij as the 
'dressed' weight of Note that G is always connected 

because H is irreducible. This completes our first task 
in identifying the appropriate graph and Laplacian for 
a given Hamiltonian, H. Next, we show how weighted 
graph can be used to give a lower bound to the spectral 
gap of the Hamiltonian. 

Cheeger constant and Cheeger inequalities: We 
now introduce the appropriate Cheeger constant associ- 
ated with a graph G. Consider a bipartition S and S 
in the vertex set V of the graph. Let the flow and the 
capacity be defined as 
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respectively. We define the Cheeger constant as 

F s 



C s <l/2 



(5) 



(6) 



where the minimisation is over all possible partitions. 
The edge boundary dS that corresponds to $ is called 
here the Cheeger cut and it identifies the bottleneck of 
the graph, as depicted in Fig. [T] 

Next we present the Cheeger inequalities jS] . Consider 
the Laplacian operator |2]) with spectral gap A(L) = 
Ai — Ao, defined on the graph G with Cheeger constant 
$. Then the Cheeger inequalities are give by 



(a) 2$ > Ai - A and (b) A x - A > 



$ 2 

2N 



(7) 




FIG. 1. An unweighted graph G with vertices V and edges 
E. The Cheeger cut, dS, is depicted that splits the vertices 
V into S and its compliment S. 



The Cheeger inequalities bound from above and below 
the spectral gap of the Laplacian and of the correspond- 
ing Hamiltonian. It is worth noting that the bounds are 
functions of the Cheeger constant $ and the 'bare' de- 
gree | Ao | which are geometric characteristic of the graph. 
Hence, if we had the means to determine the general 
shape of the graph from general properties of the Hamil- 
tonian then we could successfully estimate the energy gap 
A(ff). 

The proof of the upper bound (J7^,) is rather simple and 
it is based on a variational argument for the first excited 
state of the Laplacian [5]. By taking a vector tp that 
is orthogonal to the ground state of the Laplacian we 
guarantee that the expectation value of L with respect 
to ip will be larger or equal to the gap Ai — Ao- Such a 
vector can be defined on a bipartition of G into A and B 
so that 




(8) 



The variational parameter is the position of the boundary 
that separates A and B. The optimal value is obtained 
when the boundary is identical to the Cheeger cut, where 
the gap is upper bounded by 2$ as shown in (JT^i). Hence, 
Hamiltonians that give rise to graphs with a predominant 
bottleneck behaviour, i.e. small $, have a small gap. If 
on the other hand $ is large it does not automatically 
guarantee that the Hamiltonian has a large gap. This is 
due to the factor 1/|A | in the lower bound (|7p), where 
|Aq|, roughly speaking, corresponds to the size of the sys- 
tem. As we are interested in the behaviour of the lower 
bound in the thermodynamic limit we are faced with the 
task of improving this bound. 

Two characteristic examples: We shall now see 
with explicit examples how the Cheeger inequalities are 
applied to physical systems. Our first example concerns 
a free particle hoping on a one dimensional lattice. For 
simplicity we take the lattice to have periodic boundary 
conditions and to be of size N. The Hamiltonian is given 
by 



(9) 
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FIG. 2. The graph corresponding to a particle in one di- 
mensional lattice of N sites with hopping amplitude t subject 
to periodic boundary conditions. For this system the graph 
as well as the physical lattice configuration are the same. 
Symmetry considerations show that the Cheeger constant is 
$ = At/N. 



(a) (b) (c) (d) 




FIG. 3. The hypercube graphs Q n for (a) n = 1, (b) n — 2, (c) 
n — 3 and (d) n = 4. In all cases one can see that the optimal 
cut, dS, has S and S being lower dimensional hypercubes 

Qn-l- 



where the state |i) denotes the particle being in position 
i with the periodic condition \N + 1) = |1). The ground 
state is given by \ipo) = h= J2i=i I*) with eigenvalue 
Ao = — 2t. By Fourier transformation one can show that 
Ai = —2t cos(ir / N) . Hence, in the thermodynamic limit 
A(H) = 2t[l - cos(tt/A0] tn 2 /N 2 -> as iV -> oo. 

Let us now turn to the graph theoretic approach of this 
system. Following definitions ( 2|3|4 ) we obtain the graph 
of Fig. |2]with Cheeger constant $ = At/N. The Cheeger 
inequalities give 



8t 4< 

> Ai - A > , 



(10) 



that is in agreement with the exact result. Hence, the 
estimation of the energy gap from the Cheeger inequal- 
ities gives the gapless asymptotic behaviour of the one 
dimensional free particle when N — > oo. 

Next we consider the Hamiltonian of non-interacting 
spin-1/2 particles 



(11) 



where a x is the Pauli operator and B > 0. The ground 
state is an equal superposition of all spin states |V>o) = 
Y^i=i K) where \i) denotes a particular spin config- 
uration among the total of N — 2" in the a z eigen- 
basis. The two lowest eigenvalues are Ao = —Bn and 
Ai = —B(n — 2) giving the gap A(H) = 2B, which is 
constant as n is taken to the thermodynamic limit. 

The corresponding graph is a hypercube Q n depicted 
in Fig. [3j The weights of the edges are all B so that the 
Cheeger constant is $ = B. The Cheeger inequalities 
give 



2B > A(H) > 



B 
2n' 



(12) 



We see that, while both inequalities hold, the lower 
bound goes to zero when n — > oo. Hence, the Cheeger 
inequality gives no information if the Hamiltonian H is 
gapped or not in the thermodynamic limit. This simple 
example demonstrates that the current lower bound is 



not adequate for our purposes. In the following we intro- 
duce a generalisation of the lower bound of the Cheeger 
inequalities that can recognise if a Hamiltonian is indeed 
gapped or not in the thermodynamic limit. 

Generalised Cheeger inequality: The inadequacy 
of the lower bound (J7Jd) comes from the factor l/|Ao|, 
where |Ao| corresponds to the maximum number of edges 
originating from a vertex. In the following we show that 
it is not necessary to consider all edges of the graph G in 
the derivation of the lower bound. In fact, a generalised 
form of the lower bound can be derived that does not 
suffer from the scaling problem of the original version 
and, hence, it can be applied to problems of physical 
interest. 

Let us assume we know the optimal bipartition of a 
graph G in S and S by the Cheeger cut dS. We con- 
sider a reduced version of the graph G, named G, that 
has exactly the same vertices as G, but a reduced set of 
edges, E. The 'bare' degree of vertices in G are given 
by Q = J^j (j w ij/ a i - Hence, G has a smaller maxi- 
mum 'bare' degree than G namely, c = max^gy Cj < |Ao|. 
Consider now all possible subsets Si of the set of vertices 
of S. For each subset we define the ratio $^ = FsJCsi 
where i 7 s i is defined only in terms of the edges of the re- 
duced graph and we set $ = ming <I>i. Note that, unlike 
the minimisation in the definition of $ is performed 
only over the set S not over the whole G. Nevertheless, 
since the reduced graph has a smaller number of edges, 
the value of $ might be actually smaller than 

With the help of the reduced graph one can derive (see 
Supplement) a generalised version of the lower bound: 



Ai - A > 



$ 2 

2c' 



(13) 



This inequality shows that the energy gap of H can be 
bounded by geometric characteristics of a reduced graph 
G. Initially, note that for G = G the lower bounds (17p) 



and ( 13 ) are identical. As the reduced graph can be arbi- 



trarily chosen we aim to define G in order to maximise the 



lower bound (13 1. If a reduced graph can be found with 



a lower Cheeger bound that does not tend to zero in the 
thermodynamic limit then the Hamiltonian is gapped. 
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FIG. 4. Possible reduced graphs for (a) the general graph of 
Fig. I] and (b) for the Qi hypercube. These reduced graphs 
are chosen to optimise the lower Cheeger bound. The removed 
edges are denoted in light grey. 



The advantage of this approach is that for certain 
Hamiltonians we can easily make c finite as the size of 
the system, N, increases. This can be performed by hav- 
ing only a fixed 'bare' degree in G that does not increase 
with N. Moreover, it is desired to have \dSi\ ^ for all 
subsets Si so that the inequality ( 13 ) does not become 
trivial. To meet these requirements a general rule can 



be adopted of keeping in E all edges that belong in the 
Cheeger cut and paths of edges that connect all vertices 
in S with vertices that have edges in dS (see for example 
Fig.ga)). 

Even if at a first sight the choice of the reduced graph 
might seem random it is rather straightforward to demon- 
strate its versatility. Consider, for example, the case of 



the simple spin-1 /2 model given in (11). The energy spec- 
trum of this model is gapped in the thermodynamic limit, 
but the usual Cheeger inequalities failed to demonstrate 
this. Let us define the reduced graph such that all edges 
are removed except the ones that belong to the Cheeger 
cut (see Fig. mb)). This has c = B and $ = <f> = B 
giving the lower bound 



Ai - A > 



B 



(14) 



which is independent of the system size. Hence, the 
generalised Cheeger inequality can successfully prove the 
gapped nature of this model in the thermodynamic limit. 
Consider now the one dimensional Ising model 




FIG. 5. The upper and lower bounds of the energy gap for 
the Ising model, (a) The spectral gap (solid line) for n = 14. 
The upper and lower Cheeger bounds (dotted lines) and the 
generalised lower bound (dashed line) as a function of the 
magnetic field, B. (b) The lower Cheeger bound (dotted line) 
and the generalised lower bound (dashed line) as a function of 
the system size n for B = 2. It is clearly seen that, unlike the 
usual Cheeger bound, the generalised one successfully bounds 
the gap as the system size increases. 



be straightforwardly applied to determine theoretically 
or numerically the gapped nature of more complex sys- 
tems [IT] . 

Conclusions: In this Letter we generalised the vari- 
ational method that determines an upper bound for the 
energy gaps of Hamiltonians to the case of lower bounds. 
To upper bound the energy gap one employs the vari- 
ational ansatz for the first excited state. This method 
introduces variational parameters as degrees of freedom 
with which one can optimise the bound. Similarly, in 
the case of the lower bound choosing the reduced graph 
G provides the freedom for optimisation. In the pres- 
ence of ground state degeneracy due to a symmetry one 
can project the Hamiltonian to irreducible components 
and treat separately each of them. We envision that the 
presented method can actually reveal the behaviour of 
the gap in the thermodynamic limit of a wide variety 
of physical problems ranging from condensed matter to 
high energy physics. 
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As this is a spin- 1/2 Hamiltonian it gives rise to a hyper- 
cube graph. At B = it has a doubly degenerate ground 
state, at B = 1 it exhibits critical behaviour as n — > oo 
and it has a unique gapped ground state at B > 1. Fig. 
[5]ja) illustrates the usual and generalised Cheeger bounds 
for various values of the mag netic field B. Fig. gb) 
shows that the usual Cheeger bound scales unfavourably 
as a function of n, while the generalised one remains al- 
most constant as n increase, successfully bounding the 
energy gap from below. This general prescription can 
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SUPPLEMENT 

Recall that for a Hamiltonian H, the corresponding 
Laplacian is given by 

L = -\ \ N + D- 1 HD. 

Let %l>\ be an eigenvector of H with eigenvalue Ai, so e = 
tjjiD = (ei)^g 1 is a left eigenvector of L with eigenvalue 
Ai < Ao- Then 

eL = (Ai - A )e. (16) 
Define the subset of vertices 

V + = {i G V : e t > 0}. 
Note that J2 ieV &% = because 

0=(V>i,^o) = (V>iA£>-Vo) = (e,l)=^ ei . 

We may assume without loss of generality that CV+ = 
Sigs a i — \- Now let e be the vector defined by 

\ei/wi, i G V + ; 
I 0, otherwise. 

We assume that eo > &\ > . . . ejy-i, which implies that 
V + = {0, 1, ... , r} for some < r < N - 1. Taking the 
scalar product of ( 16 ) with e gives 



since both Ljj and are negative for j <£ V + ; 



{eL, e) = (A a - A )(e, e) = (Ai - A ) ^ tt,^. (17) 

The left hand side can be expanded as 

(eL,e)= J2etL Jt ej> ^ ^ <'./•,-'.;• (18) 



Yl ei(^T H ii) e i- 5Ie l (^^)e l , 



iev+ jev + 



(19) 

by Lji = ^H 3i and L u = -A + = - Ei/i ^7^; 
- = -S X] ? %- g i e i + XJwye?, (20) 

by ;/', , = n,//,,d,: 

... = -2^u; ij -eie J - + ^u;y(e? + e|), 



(21) 



by using j = 2 Ei<j; 



(22) 



Putting together (17) and (22 1 we get 



A(L) = Ar - A > - i<J v 



Y Ji<i w l3 {e t -e J f 



(23) 



Now we are ready to prove the following theorem. 

Theorem 1. For the the Laplacian L given in |1|) it 
holds that 



Ai - A > 



$2 

2c' 



where the quantity $ is defined so that the reduced flow 
out of S , Fs satisfies 



F s > ®C S 



for any S C V and, 



Fs = Wij, the reduced flow out of S; 

i<£S,j£V 

c = max( Wij/ Hi), the constriction. 



Proof: To prove Theorem 1 we employ the max- flow 
min-cut theorem [7]. Unlike the standard proof that is 
based on a specific vertex enumeration [5], the max- flow 
min-cut theorem provides a much more versatile frame- 
work that allows us to successfully generalise the lower 
bound. Consider the network N, based on the graph G, 
with vertex set {s,t} U X UY where s is the source, t 
is the sink, AT is a copy of V + and Y is a copy of V. 
The directed edges of this network and their capacities 
are given as follows: 
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1. For every x £ X, the directed edge (s, x) has ca- 
pacity (1 + $)7Ti where x is labelled by vertex i. 

2. For every x £ X, y £ Y, there is a directed edge 
(x, y) with capacity wiy if x is labelled by a vertex 
i, y is labelled by a vertex j and is an edge 
in G such that € E. Otherwise we set the 
capacity of to be zero. 

3. For every x £ X, y £ Y, there is a directed edge 
(x, y) with capacity 7Tj + u>u if x and y are labelled 
by the same vertex. 

4. For every y £ Y labelled by j, the directed edge 
(y, t) has capacity 7T,-. 

We claim that the value of the min-cut of this network 
is (1 + $)Cy+, where C v + = Eiev+ w i- ^° snow this, 
let .ft' denote a cut separating s and i. Let 

Xt : = {x G X : (s,x) ^ K} 
Y 1 : = {y£Y :(y,t)£K}. 

Define the following sets of edges, 



Recall (23) 



Ex 

E{X l ,X 1 ) 



= {(,s,x) : x G X}, 
= {(s,x) : x G X x }, 

= {(x, y) : x £ Xi C X and?/ G Xi C F}, 



The total capacity of the cut K is at least the sum of 
capacities of the edges, 

• Edges E x \ E Xl , 

• The edges E(X 1 ,X l ), 

• The edges E{X U Y). 

Therefore, the total capacity of the cut is at least equal 
to 

(l + *)C X - Xl +C Xl +F Xl , 
> (1 + $)C X -x 1 + C' Xl + $C Xl , by F s > 
= (1 + $)C X . 

Since there is a cut of size (1 + <fr)Cy+ , we have proved 
that the min-cut is equal to (l + $)Cy+ . By the max-flow 
min-cut theorem [7], there exists a flow function hij for 
all directed edges in the network so that is bounded 
above by the capacity of (i, j). Also, for each fixed x £ X 
and y £ Y, we have 



J2 h *J = (1 + $)tTz> 
jeY 



(24) 
(25) 



Ai - A > 



E,v, - 



and note that 



Ax - A > 



EiGV+^e? E *<j w ij (e l + e J ) 2 



At this point we have introduced the reduced graph G 
through the summation over the edge set E. This has 
been possible as we are allowed to insert any number in 
the right hand side of the above inequality as long as it 
is smaller or equal to one. Exactly this mathematical 
freedom allowed us to generalise the Cheeger inequality. 
In the denominator 



E w i j{ei + e j f<2 w ij(^ 
{hj}eE 



2 1 e 2 ) (26) 



w^ + e 2 ) (27) 

2 E^ 2 E w v ( 28 ) 

i€V j 



< 



2J2 

iev 

2C X TTi 

i£V + 



In the numerator 



i<j 
E 



i£V+ 



> 



( E ^) 2 * 2 



(29) 
(30) 

(31) 



E E 

h 2 

>Y /WlJ (e l -e J ) 2 Y / ^(e l + B ] ) 2 (32) 

E E 

> (^2 hij(e 2 - e 2 )) 2 , by Cauchy-Schwarz; (33) 

i<j 
E 

= (J2^(e 2 ~e 2 )) 2 , (34) 



= ( E M^ 2 ~e 2 )) 2 , by using £ = (35) 

{i,j}eE i<j 

= (E^ 2 ( E ^- E M) 2 (36) 

iev j j 

{i,j}£E {i,j}eE 
>(]T ^((l + ^-TT,)) 2 



(37) 
(38) 



Putting ( 38 1 and ( 30 ) together we get 



as required. 



